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INTRODUCTION  AND  RELEVANCE  OF  THE  PROBLEM 

The  recently  increased  interest  in  space  for  civilian  and  defense 
purposes  has  resulted  in  new  concepts  for  the  hypersonic  vehicles.  The 
Aeroassisted  and  Aeromaneuvering  Orbital  Transfer  Vehicles  (Refs.  1-2)  (AOTVs 
and  AMOTVs),  the  Transatmospheric  Vehicles  (TAVs)  and/or  Orbit-On-Demand  (OOD) 
vehicles  (Ref.  3),  the  Hypersonic  Cruise  and  Manuever  Airplanes  (HCMA)  (Ref. 
4),  and  the  present  (Ref.  5)  and  next  generation  space  shuttles,  are  a  few  of 
the  examples  of  these  vehicles  which  have  created  a  resurgence  of  interest  in 
hypersonic  aerothermodynamics . 

Since  most  of  the  future  hypersonic  vehicles  will  be  operating  in  the 
upper  atmosphere,  the  'low  density  effects'  will  play  a  major  role  in 
establishing  the  lift,  drag,  moments,  and  aerodynamic  heating  of  a  hypersonic 
vehicle.  In  addition,  the  new  vehicular  concepts  will  require  more  accurate 
knowledge  of  hypersonic  aerothermodynamics  under  the  low  density  conditions. 
The  so-called  'shock  and  surface  slip  effects,'  multicomponent  nonequilibrium 
gas  chemistry,  proper  modeling  of  gas-surface  interactions,  reaction-rate 
coefficients  for  high  temperature  air  as  well  as  a  correct  description  of 
thermodynamic  and  transport  properties  would  be  the  pacing  research  activities 
for  an  accurate  prediction  of  the  aerothermal  environment  of  the  new 
generation  of  hypersonic  vehicles. 

The  problems  of  hypersonic  aerodynamics  without  the  low  density  effects 
have  been  analyzed  in  detail  (Refs.  6-17)  both  for  perfect  gas  as  well  as  for 
air  with  real  gas  effects.  Not  much  attention  seems  to  have  been  given  to  the 
problems  encountered  with  low-density  aerothermodynamics.  Reference  18 
analyzed  the  flow  past  a  flat  plate  at  low  densities  with  body  slip  only. 
Davis  (Ref.  19)  included  body  and  shock-slip  in  the  viscous  shock-layer 
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analysis  of  perfect-gas  flow  around  a  hyperboloid.  The  low-density  shock, 
thickening  effects  with  the  integration  of  full  Naviei-Scokes  equations  in  the 
stagnation  region  were  recently  obtained  in  Ref.  20. 

The  surface  slip  conditions  obtained  by  Davis  (Ref.  19)  contained  some 
errors  as  explained  in  Ref.  21.  The  shock-slip  boundary  condition  of  Ref.  19 
also  did  not  account  for  the  derivatives  of  the  shock  quantities  in  the  shock- 
oriented  coordinate  system.  This  introduces  significant  errors  in  analyzing 
flows  past  slender  bodies  as  compared  to  the  wide-angle  bodies.  Tree,  at ■ al . 
(Ref.  22),  analyzed  the  hypersonic  ionizing  viscous  shock-layer  flow  past 
axially  symmetric  bodies  at  low  densities.  Their  analysis  was  for  an 
equilibrium  Catalytic  wall  and  contained  the  errors  of  Ref.  19  in  the 
specification  of  body-  and  shock-slip  boundary  conditions.  Lwaminathan,  et 
al.  (Ref.  23)  and  Song,  et  al.  (Ref.  24)  recently  included  the  surface-  and 
shock-slip  effects  for  the  three-dimensional  flows.  Their  Bhock-slip 
formulation  contained  the  error  mentioned  earlier,  whereas  their  surface-slip 
model  was  taken  from  Ref.  25,  which  also  contained  errors  as  explained  in  Ref. 
21.  Further,  even  without  the  shock-  and  surface-slip  boundary  conditions, 
the  viscous  shock-layer  (VSL)  three-dimensional  code  of  Refs.  23  and  24  gives 
flowfield  quantities  (such  as  boundary  layer  thickness,  flowfield  profiles  of 
velocity,  temperature,  etc.),  which  are  physically  not  very  appealing  (Ref. 
26).  The  surface  heat  transfer  and  skin-friction  coefficient  predicted  by  the 
code  developed  by  Clark  Lewis  and  his  associates  (Refs.  23  and  24)  is  impacted 
considerably  by  these  poor  flowfield  quantities  for  a  turbulent  flow  as  shown 
in  Ref.  26.  The  inaccurate  flowfield  profiles  have  also  resulted  in  poor 
comparison  of  predicted  surface  heat  transfer  with  the  shuttle  flight  data 
(Ref.  27).  One  of  the  major  drawback  with  Lewis'  code  is  that  it  smooths  the 


shock-stand-off  distance  after  each  global  flowfield  solution,  resulting  in 
poor  conservation  of  mass  and  momentum  in  the  flowfield.  Instead,  in  the  code 
developed  here,  the  streamwise  derivatives  of  the  shock  stand-off  distance  are 
smoothed.  This  results  in  superior  conservation  of  mass  and  momentum  in  the 
flowfield . 

An  additional  premium  is  placed  on  the  theoretical  analysis  by  the  fact 
the  most  of  the  ground-based  experimental  facilities  fail  to  simulate  the  low- 
density  highly  energetic  flowfields.  In  particular,  the  interaction  between  a 
high  enthalpy  gas  flow  and  a  relatively  cold  surface  for  the  low-density 
hypersonic  re-entry  conditions  is  poorly  simulated  (at  best)  with  the  existing 
experimental  facilities.  These  limitations  in  the  simulation  by  ground-based 
facilities  resulted  in  the  design  of  the  Galileo  heat  shield  (Refs.  14-17,  28) 
completely  by  means  of  detailed  flowfield  calculations.  This  is  considered  a 
'first'  in  modern  hypersonics  where  no  experimental  data  in  the  final  design 
of  heat  shield  were  used  since  they  were  not  available.  Thus,  there  is  a  need 
to  develop  an  accurate  viscous  shock-layer  code  which  can  give  physically 
appealing  results  over  a  wide  range  of  flow  conditions. 

The  viscous  shock-layer  (VSL)  analysis  (Ref.  11)  is  one  of  the  most 
extensively  used  continuum  computational  techniques  for  analyzing  the 
flowfields  past  two-dimensional  and  axisymmetric  bodies  of  revolution  at 
hypersonic  velocities.  The  method  of  reference  11  has  been  used  both  for 
perfect  gas  as  well  as  for  air  with  finite-rate  chemistry  to  allow  for  the 
real  gas  effects.  Most  of  the  results  obtained  so  far  with  the  VSL  method  of 
Ref.  11  do  not  allow  for  the  low-density  effects  and  the  shock  is  treated  as  a 
thin  discontinuity.  The  main  advantages  of  the  VSL  method  are  that  even  with 
multicomponent  gas  chemistry,  it  can  give  results  around  the  complete  body  in 
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a  very  small  computational  time  and  also  the  computer  storage  requirements  are 
relatively  small.  The  major  deficiencies  of  the  VSL  code  of  Ref.  11  are  that 
it  does  not  account  for  the  low-density  effects  as  well  as  it  cails  to 
converge  for  slender  bodies  (such  as  sphere  ^ones  with  body  half-angles  of 
less  than  25°). 

Reference  20  employs  the  complete  Navier-Stok.es  equations  where  no 
assumption  is  made  about  the  shock  thickness  and  the  shock  is  captured  as  part 
of  the  solution.  The  main  disadvantage  with  this  approach  is  large  storage 
requirement  because  the  numerical  integration  is  carried  out  to  the  free- 
stream.  The  computational  times  will  also  vary  large  if  solutions  are  desired 
around  a  complete  body.  These  solutions,  however,  can  serve  the  basis  for 
changes  and/or  improvements  required  in  the  VSL  technique  under  low  density 
conditions. 

In  comparison  to  the  noncontinuum  methods  like  the  Direct  Simulation 
Monte  Carlo  (Ref.  29)  (DSMC),  etc.,  the  continuum  methods  (such  as  VSL  and 
Navier-Stokes)  are  fast  and  give  practical  results  in  a  fraction  of  the 
computational  time.  The  recent  work  (Refs.  5,  20,  21)  shows  the  accuracy  of 
continuum  prediction  techniques  if  proper  accounting  is  made  of  the  low- 
density  effects. 

As  mentioned  earlier,  the  experimental  data  for  the  low-density  high 
energy  flows  are  almost  non-existent  in  the  literature.  For  low  energy  flows 
and  for  shapes  with  spherically-blunted  nose-tips  (such  as  hyperboloids  and 
sphere  cones),  Refs.  30  and  31  have  provided  good  data.  Little's  data  (Ref. 
30)  are  limited  to  the  measurements  of  pressure,  drag,  and  skin  friction, 
whereas  Boylon's  data  (Ref.  31)  are  for  the  surface  heat  transfer.  These 
measurements  were  made  during  the  nineteen-sixties  and  early  nineteen- 
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seventies.  Thus,  there  is  a  considerable  need  to  have  the  experimental,  in 
particular,  surface  heat  transfer  data  under  the  low  density  flow  conditions. 
With  the  superior  data  acquisition  and  reduction  techniques  available  now,  it 
should  be  possible  to  obtain  the  new  data  with  fairly  high  degree  of  accuracy. 

PRESENT  WORK 

With  the  deficiencies  and  inadequacies  of  the  continuum  technqiues 
(outlined  in  the  introduction)  in  mind,  a  viscous  shock-layer  (VSL)  code  has 
been  developed  ror  the  low  density  flight  conditions  for  wide-angle  as  well  as 
long  slender  bodies.  Under  the  low  Reynolds  number  (or  low-density)  flow 
conditions,  the  viscous  effects  influence  almost  the  entire  shock-layer  and 
the  shock  itself  is  considerably  thick  as  compared  to  the  high  Reynolds  number 
(or  high-density)  case.  This  requires  the  use  of  more  comprehensive  set  of 
governing  equations  as  compared  to  the  classical  boundary-layer  equations. 
The  complete  Navier-Stokes  equations  are  considered  (Ref.  20)  appropriate  for 
the  low  Reynolds  number  applications.  But  computer  storage  and  computational 
time-wise  they  are  very  expensive  to  solve  for  flows  around  the  long  bodies. 
The  viscous  shock-layer  (VSL)  equations  (Refs.  11,  16,  19)  represent  an 
intermediate  level  of  approximation  between  the  boundary-layer  and  Navier- 
Stokes  equations.  In  obtaining  the  viscous  shock-layer  equations  from  the 
Navier-Stokes  equations,  terms  are  kept  up  to  second-order  in  the  inverse 
square  root  of  the  Reynolds  number  from  both  a  viscous  and  inviscid  viewpoint. 
This  results  in  simplified  governing  equations,  which  are  uniformly  valid  to 
moderately  low  Reynolds  numbers.  To  the  order  of  approximation  involved  in 
obtaining  the  VSL  equations,  the  body  surface  conditions  are  given  by  slip  and 
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temperature  jump  conditions  (Ref.  21)  while  a  set  of  shock-slip  conditions* 
are  used  to  determine  conditions  behind  the  shock. 


ANALYSIS 

Flow  Governing  Equations 

The  conservation  equations  employed  in  the  present  analysis  are  the 
steady  perfect-gas  viscous  shock-layer  equations  (Refs.  11  and  19)  for  an 
axisymmetric  or  two-dimensional  body  at  zero  angle  of  attack  (Fig.  1).  These 
equations  in  the  orthogonal,  body-oriented  transformed  coordinates  and 
nondimensional  form  are  provided  here.  The  second-order  partial  differential 
equations  are  written  in  the  following  form: 


g2w  (3zg/Sn2)  +  ajOg/an)  aw  «2 

an2  Og/3n)2  3n  Og/9n)2 

a3  ,  °4  aw 

+  - IT-o  + - -~2  TF  “  °' 

(3g/3n)  (3g/3n>  ^ 


The  quantity  W  represents  u  in  the  s-momentum  equation  and  H  in  the 
enthalpy  energy  equation.  The  coefficients  to  are  written  as  follows: 
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These  are  similar  to  those  of  Ref.  19  except  for  the  corrections  to 
tangential  velocity  and  temperature  slip  expressions. 
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Energy,  W  *  H: 
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The  remaining  first-order  equations  are  written  as  follows: 


Global  Continuity: 
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sh 


(4) 


n-momentum: 


pu 


(i  +  n  ngh  k) 


3v 


n8h  n  PU 


n  ,  (1  +  n  n  , 
sh  sh 


__  _3v  _3g  ,  pv  3v  3g 

.  _  —  ~  n  ,  3rt  — 

ic)  3n  3n  sh  3n 


2 

pu  K 


(i  +  r,  nsh  K) 


1 _ _3p  il 

nsh  3n  3n 


0 


state : 

p  -  pT(y-l)/y. 


(5) 


(6) 


In  Eqs.  (2)  through  (5),  prime  denotes  differentiation  with  respect  to  £. 
Further,  Eqs.  (1)  through  (5)  contain  the  following  independent  transformation 
(Ref.  20): 


h  *  g(^ — )  -  g(n). 
sh 


(7a) 


The  stretching  function  g(n)  is  given  by 


9 


g(n) 


1 


a  + 


(1-a) 


t„  ji  -.iw  i  i 

(  g  +  n(  2  ot+1 )  -  1 


and  Its  first  and  second  derivatives  are 


(7b) 


dg  _  (l-a)(2orH)  j _ 1 _  +  _ 1 _ ) 

dn  '"(IS) 1  [I  -  n(2ct+l)  +  1]  [B  +  n(2orfl)  -  1]  I 


d2g  _  ( l-q)(2cH-l )2  |  1 _ 

dq2  lp»  ‘  nCTorH)  +  l]2 


_ 1 _ 

[?  +  "q(  2  ori-L  )  -  1]' 


(7c) 


(7d) 


Equation  (7b)  permits  the  mesh  to  be  refined  either  near  the  body  ("a  *  0)  or 
refined  equally  near  both  the  body  and  bow  shock  (ci  -  1/2)  when  6hock 

becomes  thick  under  the  low  density  flight  conditions.  Parameter  ~B  controls 
the  amount  of  refinement  with  values  near  1  giving  the  largest  amount  of 
stretching.  Equation  (7b)  may  be  inverted  to  obtain  the  physical  coordinate 
n  from  the  transformed  coordinate  q: 


(7e) 


The  transformation  of  Eq.  (7a)  keeps  the  body  at  q  -  0  and  the  shock  at 
q  *  1  with  uniform  mesh  in  the  computational  coordinate  q. 
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For  ease  In  numerical  computations,  Eqs.  (1)  through  (6)  are  further 
transformed  by  normalizing  most  of  the  variables  with  their  local  shock 
values.  The  normal  velocity  v  is  not  normalized  because  this  quantity 
changes  direction  at  the  shock  and  can  be  near  zero.  Therefore,  the 
normalized  v-profiles  may  not  be  well  behaved. 

Boundary  Conditions 

The  following  wall  and  shock-boundary  conditions  are  employed. 


(8a) 


(8b) 


(8c) 


P 


sh 


sin  a 


(9a) 
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(il)  s-momentum 


e2psh  |  [c°s(.a-e)  +  n  — ; ■-  sin(a-e)]  -  — 

n  3n  3n 


H 


3u 


) 


sin(a-8)jg^  +  ush  s*n  a  "  s*n  a  cos  a 


(9b) 


(iii)  n-momentum 


Psh  -  Poo  +  Sin  a(Sin  a  +  Vsh) 


(9c) 


(Iv)  Energy 


.2  (V  ■>  ir _ ✓  «•>  ,  —  8h 


'fe) «k  [cos(a-e)  +  n 


.  ,  -,-,1  3g  3T 

Pr  Jsh  )L— v-  jg—  sln(a-e)]^  —  — 

dr]  3ri 


3T  sin  a  ~  N2 

-  ^  sln(o-e)  j  gh  +  Tgh  sin  a  -  -j-  -  cos  a) 


sin  a 


4y 


sin  a  + 


(Y+ir 


_  »(T-1) 

(y+D2J 


M 


2  4  2 

( y+1 )  M  sin  a 
00 


(v)  Equation  of  State 


(9d) 


“sh  •  1  psh/(’r-1)T8h- 


(9e) 


In  Eqs.  (9a)  through  (9e),  u  and  v  are  the  components  of  velocity  tangent 
and  normal  to  the  shock  interface,  respectively,  and  are  related  to  the 
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components  in  the  body-orientcu  coordinat  co  ab  • 

%h  ’  Ush  C0S<°-6)  +  V8h  8in<a-0>  <9f) 
vgh  *  "  ugh  sin(a-e)  +  vgh  cos(a-0).  (9g) 

There  are  errors  in  the  expressions  for  surface  pressure  and  temperature 
slip  (i.e.,  Eqs.  (2.7c)  and  ( 2 . 7 d ) )  of  Ref.  19.  The  pressure  and  temperature 
slip  values  are  underpredicted  there  by  the  ratio  of  K/ p  as  explained  in 
Ref.  21.  For  a  perfect  gas  this  ratio  is  the  same  as  1/Pr.  Similarly,  Eqs. 
(2.8d)  and  (2.8e)  of  Ref.  19  for  the  tangential  velocity  and  temperature  slip 
at  the  shock  contain  normal  derivatives  of  the  shock  quantities  in  the  body- 
oriented  coordinate  system  in  place  of  the  shock-oriented  coordinates  as  given 
by  the  underlined  terms  in  Eqs.  (9b)  and  (9d). 

METHOD  OF  SOLUTION 

The  method  of  solution  employed  here  is  an  implicit  finite-difference 
method  of  Refs.  11  and  19.  However,  this  method  is  implemented  here 
differently  from  Refs.  11  and  19  because  convergence  problems  were  encountered 
in  these  references  with  the  increase  in  shock-layer  thickness. 

Let  the  subscript  m  denote  the  station  measured  along  the  body  surface 
and  n  denote  the  station  measured  normal  to  the  body  surface.  The  second- 
order  equations  (1)  through  (3)  are  replaced  with  central  differences  taken  in 
the  n  direction  and  two-point  backward  differences  in  the  £  direction  at 
the  point  (m,n).  Further,  based  on  the  hypersonic  small  disturbance  theory 
for  a  slender  body,  the  first-order  global  continuity  Eq.  (4)  and  n-momentum 
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Eq.  (5)  are  solved  for  the  pressure,  p,  and  normal  velocity,  v,  in  a  coupled 
way.  The  density  in  these  equations  is  eliminated  through  the  use  of  equation 

of  state,  Eq.  (6).  Each  of  the  Eqs.  (4)  and  (5)  are  expressed  in  finite- 

difference  form  at  points  (m,  rrt-  V2  )  and  (m,  n-  V2  )  using  a  box  scheme 
discussed  by  Richtmyer  (Ref.  32).  In  these  two  coupled  equations,  p  and  v 

are  eliminated  alternatively  to  solve  for  these  two  variables.  The  finite- 

difference  forms  of  Eqs.  (1)  through  (5)  can  easily  be  solved  by  using  the 
Thomas  Algorithm. 

The  solution  is  started  at  the  stagnation  point  where  the  various 
flowfield  quantities  are  expanded  (Refs.  11  and  20)  in  terms  of  the  distance, 
£,  along  the  body  surface.  These  series  expansions  reduce  the  partial 
differential  equations  (1)  through  (5)  to  ordinary  differential  equations  in 
terms  of  n.  At  a  body  location  m,  other  than  the  stagnation  point,  a  two- 

point  backward  difference  is  used  for  the  derivative  with  respect  to  £  at 
the  point  (m,n)  as  mentioned  earlier.  This  again  gives  ordinary 

differential  equations  at  location  m  in  terms  of  n  for  equations  (1) 
through  (5).  The  finite-difference  form  of  these  ordinary  differential 
equations  (obtained  through  the  central  differences)  can  be  solved  by  using 
the  Thomas  Algorithm  as  mentioned  earlier.  Figures  2a,  2b,  and  2c  present  the 
flow  diagrams  for  obtaining  solutions  with  body  slip,  shock-slip,  and  body  and 


shock-slip,  respectively. 
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DISCUSSION  OF  RESULTS 

Some  of  the  obtained  results  by  the  present  method  are  given  in  Figures  3 
through  13. 

Reference  19  carried  out  some  of  the  early  detailed  calculations  for  the 
low  density  flow  conditions.  In  Refs.  11  and  19,  the  governing  equations  are 
solved  in  a  successive  manner  and  an  initial  shock  shape  is  required  as  an 

input  so  as  to  make  the  viscous  shock-layer  equations  (which  are  weakly 

elliptic  in  the  tangential  direction)  parabolic.  This  allows  the  use  pf  a 
numerical  scheme  with  space  marching  in  the  streamwise  direction.  In  the 

approach  of  Refs.  11  and  19,  the  initial  shock  shape  is  obtained  from  the 
thin-shock-layer  approximation,  whereas  Refs.  23  and  24  have  obtained  it  form 
an  inviscid  solution.  Both  of  these  approaches  are  unworkable  under  the  low- 
density  conditions,  especially  for  a  slender  body.  In  the  present  method  of 
solution,  the  initial  shock  shape  is  obtained  from  the  corresponding  cone 

angle  solution  (Ref.  26)  for  a  high  Reynolds  number  flow  case.  For  a  very  low 
Reynolds  number,  the  initial  shock  shape  may  be  obtained  by  reducing  the 
Reynolds  number  sequentially  by  a  factor  of  10  or  so.  This  approach  gives 
very  good  estimate  of  the  initial  shock  shape  and  results  in  consistent 
convergence  of  the  solutions.  An  alternate  method,  which  avoids  prescription 
of  the  initial  shock  shape,  has  recently  (Ref.  33)  been  suggested  for  the  high 
density  flows.  However,  this  method  is  not  suitable  for  long  slender  bodies. 
Moreover,  for  the  faster  convergence  and  stability,  a  good  initial  shock  shape 
is  desirable  for  this  method  also.  In  the  present  approach,  the  continuity 
and  normal  momentum  equations  are  solved  simultaneously  as  a  coupled  set  to 
overcome  the  stability  problems  encountered  in  Refs.  11  and  19  for  slender 
bodies.  The  tangential  momentum  and  energy  equations  are  still  solved  in  a 
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successive  manner.  As  mentioned  earlier.  Refs.  11  and  19  solve  all  the 
governing  equations  successively. 

Detailed  comparisons  of  the  predictions  with  the  experimental  data  of 
Little  (Ref.  30)  were  made  in  Ref.  19.  These  data  are  still  considered  quite 
good.  Unfortunately,  Little's  data  are  limited  to  the  measurements  of 
pressure,  drag,  and  skin  friction  only.  There  are  no  surface  heat  transfer 
data  provided.  Only  good  data  for  the  surface  heat  transfer  are  obtained  in 
Ref.  31.  Thus,  there  is  a  considerable  need  to  have  the  surface  heat  transfer 
data  under  the  low  density  flow  conditions. 

Figures  3  and  4  give  comparisons  for  the  Stanton  number  and  skin-friction 
coefficient,  respectively,  obtained  by  the  present  method  and  by  that  of  Ref. 
19.  The  calculations  are  carried  out  for  the  stagnation  point  only  for 
different  values  of  the  Reynolds  number  parameter  s,  which  i6  a  measure  of 
the  degree  of  rarefaction.  Larger  values  of  e  imply  increased  slip  (or  low 
density)  effects.  The  two  predictions  have  similar  trend.  However, 
significant  differences  are  noticeable  for  larger  values  of  e.  The  present 
calculations  give  better  agreement  with  the  experimental  data  (as  will  be 
shown  in  Figures  5  through  9)  and  are,  therefore,  considered  accurate.  The 
discrepancies  between  the  present  predictions  and  those  of  Ref.  19  may  be 
explained  by  the  following:  First,  the  transformation  employed  (in  the 

present  method)  to  cluster  the  computational  points  near  the  body  and  shock 
gives  more  accurate  slip  values.  Second,  there  are  some  errors  (see  Ref.  21) 
contained  in  Ref.  19  in  the  surface  pressure  and  temperature  slip  expressions 
as  well  as  in  the  equations  giving  slip  values  of  tangential  velocity  and 
temperature  at  the  shock  as  discussed  earlier.  The  latter  error  may  not 
affect  the  values  at  the  stagnation  point. 
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Figures  5  through  8  give  comparisons  between  the  predictions  and  the 
experimental  data  of  Ref.  30.  Clearly,  the  drag  coefficient  (with  shock  and 
body  slip)  predicted  by  the  present  method  is  in  much  better  agreement  with 
the  experimental  data  in  comparison  to  the  prediction  of  Ref.  19,  over  a  wide 
range  of  the  freestream  Reynolds  numbers,  Re,“. 

The  calculations  contained  in  Figures  2  through  8  are  for  hyperboloid 
shaped  bodies  (see  Fig.  1).  Figure  9  gives  comparison  of  heat  transfer 
distribution  for  a  ten-degree  sphere  cone.  The  present  predictions  with  shock 
and  body  slip  are  in  good  agreement  with  the  data  of  Ref.  31  except  for  the 
stagnation  point.  The  heat  transfer  rates  at  the  stagnation  point  were 
determined  to  be  biased  upward  (Ref.  34)  due  to  particle  impact  caused  by  the 
arc  heater. 

Figures  10  and  11  show  comparisons  of  predicted  skin-friction  coefficient 
and  Stanton  number  distributions.  The  comparison  between  the  present  values 
and  those  of  Ref.  33  is  quite  good.  However,  the  method  of  Ref.  33  is  fully 
coupled  as  compared  to  the  coupling  between  the  normal  momentum  and  continuity 
equations  only  in  the  present  method.  The  full  coupling  requires  solving 
a  5x5  matrix  at  every  point  in  the  flowfield  for  a  perfect  gas.  The 
complexity  and  stability  problems  (Ref.  33)  in  a  fully  coupled  solution  will 
increase  in  analyzing  a  seven-species  high-temperature  air  flow.  Also,  the 
computational  times  will  be  considerably  large  for  long  slender  bodies  by  this 
method.  Present  approach  is  more  appealing  for  such  flow  conditions.  Figures 
10  and  11  also  give  results  with  and  without  slip  for  a  variable  grid  near  the 
shock  and  body  surface.  Clearly,  the  computational  grid-size  as  well  as  the 
slip  effects  are  important  in  this  case. 
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Finally,  Figures  12  through  14  provide  a  severe  test  for  the  present 
computational  method  and  the  surface  and  shock-slip  boundary  conditions.  The 
flow  analyzed  in  these  figures  is  a  high  Mach  number  (M  -  20)  flow  over  a 
highly-cooled  (Tw  *  540°R) ,  one-inch  nose  radius,  10°  sphere  cone.  The 
surface  drag,  Stanton  number,  and  surface  pressure  distributions  are  predicted 
for  very  high  to  low  altitude  flight  conditions.  As  the  body  moves  from  low- 
density  to  high-density  flight  regime,  the  slip  effects  are  gradually  reduced 
and  become  almost  negligible  at  a  height  of  50,000  ft.  However,  they  persist 
to  a  much  lower  altitude  (namely,  50,000  ft.)  due  to  the  highly-cooled  surface 
condition  for  this  high  Mach  number  flow  case.  Also,  the  slip  effects  persist 
in  surface  drag  and  pressure  through  a  wider  altitude  range  as  compared  to  the 
Stanton  number  or  surface  heat  transfer. 

CONCLUSIONS 

Results  have  been  obtained  for  the  surface  pressure,  drag,  heat  transfer 
as  well  as  skin-friction  coefficients  for  hyperboloids  and  sphere-cone  shaped 
slender  and  wide-angled  bodies  under  varying  degrees  of  low-density  flow 
conditions.  The  flow  cases  analyzed  include  highly-cooled  surfaces  in  very 
high  Mach  number  flows.  The  present  method  is  found  to  predict  various 
flowfield  quantities  quite  accurately  when  compared  with  the  experimental 
data.  A  logical  extension  of  this  work  will  be  to  analyze  the  multicomponent 
air  flow  under  the  low-density  conditions.  The  present  results  are  limited  to 
the  flow  of  a  perfect  gas  only. 


« 
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Figure  1.  Coordinate  Syst 


INITIAL  GUESS  FOR  ALL  SHOCK  QUANTITIES 
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Figure  2c.  Shock  and  Body  Slip. 


Figure  3.  Stanton  Number  in  the  Stagnation  Point  as  a  Function  of 
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Friction  in  the  Stagnation  Point  as  a  Function  of 


Figure  6.  Drag  Coefficient  for  a  10°  Hyperboloid. 
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Figure  9.  Heat  Transfer  Distributions  Over  a  10°  Sphere-Cone. 
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Figure  11.  Stanjton  Number  Distribution  ,for  a  22.5  Hyperboloid 
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Wall  [ res  sure  as  a  Function  of  c  for  a  1C'  Sohere-Cone 


